#include "Element.h"

int main()
{
    vector<double> l_b = {-1,-1};
    vector<double> r_t = {1,1};
    size_t nx = 1<<7;
    size_t ny = 1<<7;
    int acc = 4;
    // auto f = [](double x, double y)-> double { return 2*cos(x-y);};
    // auto u = [](double x, double y)-> double { return cos(x-y);};
    auto f = [](double x, double y)-> double { return 1;};
    auto u = [](double x, double y)-> double { return -(x*x + y*y)/4;};
    Q1_Rectangel_fem_solution solution(l_b,r_t,nx,ny,acc,f,u);
    solution.solve();
    // std::cout << solution.get_Mat() <<std::endl;
    // std::cout <<solution.get_Rhs() << std::endl;
    return 0;
}